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Abstract 

We consider a scalar-metric gauge theory of gravity with independent 
metric, connection and dilaton. The role of the dilaton is to provide 
the scale of all masses, via its vacuum expectation value. In this theory, 
we study the renormalization group flow of the dilaton potential, taking 
into account threshold effects at the Planck scale. Due to the running 
of the VEV of the dilaton all particles that would naively seem to have 
masses larger than Planck's mass, may actually not propagate. This 
could solve the problem of unitarity in these theories. 
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1. Introduction 



The general ideas that go under the name of "effective field theory" are playing an increas- 
ingly important role in elementary particle physics. The variety of physical phenomena is 
divided into energy ranges, whose boundaries usually coincide with the masses of various 
particles. In each range one has an effective field theory. Very often the theories describ- 
ing neighbouring energy ranges are of the same type, the only difference being that one 
particle which is present in the higher energy range has been "integrated out" in the lower 
energy range. In other cases, however, the description of the physics in one energy range 
is quite different from that in the next energy range. It is clearly desirable to relate all 
these descriptions, but this has not always been possible so far. 

It is widely believed that Einstein's General Relativity is also an effective theory, 
valid for energies lower than Planck's energy [1,2]. This does not mean that it can only be 
treated classically; it only means that any quantum calculation in Einstein's theory will 
have a natural cutoff at the Planck scale. Prom this point of view, the nonrenormalizability 
of Einstein's theory is not a problem. 

Just below Planck's energy corrections due to higher derivative terms are expected to 
appear. The most general action with at most four derivatives is 



where R^u^ a denotes the curvature tensor constructed with the metric g^j^. As emphasized 
in [1] the gravitational action (1.1) can be treated as a quantum field theory using the 
rules of "chiral perturbation theory", which were devised in the context of a theory of 
mesons [3]. If regarded as "fundamental", the theory (1.1) is renormalizable [4,5], but 
another problem appears: the terms quadratic in curvature contain higher derivatives of 
the metric and therefore violate unitarity. Again, this is not a problem in the "effective 
theory" picture: the ghosts have masses of the order of Planck's mass, and therefore will 
never be excited if one remains at lower energies. At energies exceeding the Planck energy, 
some "new physics" is expected to appear. 

As mentioned above, there is no compelling reason to believe that the theory describing 
the physics above the Planck energy will be of the same type of the theory (1.1): it may be 
described by a completely different set of degrees of freedom. Nevertheless, in this paper we 
will conservatively assume that the "new physics" can still be described by a field theory. 
Motivated by the success of the gauge principle in explaining all other known forces of 
nature we consider a gauge theory of gravity, with independent metric and connection. We 
will assume that the action is at most quadratic in curvature and torsion. These theories 
have a long history [6]; from a particle physicist's point of view their most attractive 
feature is perhaps that they present a gravitational analogue of the Higgs phenomenon 
[7]: the vierbein behaves like a Higgs field and when it acquires a nonvanishing, constant, 
vacuum expectation value, its kinetic term (torsion squared) becomes a mass term for the 
connection. The mass is obviously of the order of Planck's mass, so below Planck's energy 
the connection degrees of freedom cannot be excited. Yet the connection itself does not 
vanish: it becomes the Levi-Civita connection, whose components in a coordinate frame 
are the Christoffel symbols. In this way the theory we will consider gives rise naturally to 



I 



d'^x ^/detg A + 



IQttG 



1 



R + aiR^^p^R''''P'' + a2R^uR^'' + a:^R'' , (1.1) 
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the action (1.1) as an effective theory below the Planck scale. The descriptions of physics 
above and below Planck's energy are easy to relate in this case. The occurrence of the Higgs 
phenomenon is often related to unification, and indeed the theory can be generalized so as 
to unify gravity with all other interactions (assuming that these are already put together 
in a GUT). The vierbein then appears as the order parameter that breaks the symmetry 
between gravity and the rest [7]. This is a true unification, in the sense in which the word 
is used in particle physics, and does not require higher dimensional spacetimes. 

The problem with these gauge theories of gravity is that there seems to be no one of 
them which is at the same time unitary and renormalizable. This is the same dilemma that 
one faces in (1.1): without the higher derivative terms the theory is not renormalizable, and 
with them it is not unitary. However, this time the problem cannot be avoided because 
we are not aware of any higher mass scale which might act as a cutoff and allow these 
theories to be regarded as effective field theories. Any theory that purports to describe 
physics above the Planck scale had better be consistent on its own. 

Our point of view is that we do not know enough on these theories to draw any 
conclusion. The reason is, obviously, that they are very complicated. In the general case 
we only have a tree level analysis of their propagators [8]. A more rigorous study of the 
spectrum would be desirable, but for the moment has been done only for Einstein's theory 
[9]. It is quite possible that the quantum propagators are very different from what they 
seem to be at the tree level. For example, some degrees of freedom may be confined. To 
the best of our knowledge, no one has performed any quantum calculation so far. 

Since the ghosts and/or tachyons would have masses of the order of Planck's mass, the 
issue of unitarity could be resolved in a rather drastic way if particles with such masses 
were forbidden to propagate. In this paper we elaborate on our earlier proposal for a 
mechanism in which this could happen. There are two main ingredients in this proposal: 
the dilaton, and the renormalization group. The dilaton is a scalar field, coupled to the 
metric, connection and to other matter fields in such a way that all masses are equal to 
its vacuum expectation value (VEV), multiplied by some dimensionless coupling constant. 
This field is closely associated to the (quantum mechanical) breaking of Weyl invariance in 
a manner which has been discussed in detail in [10,11] and is reviewed in Appendix A. In 
the present work, we compute the renormalization group fiow of the dilaton potential. In 
particular, we are interested in the running of the VEV of the dilaton, since this gives the 
dominant contribution to the running of the masses. (Dimensionless coupling constants 
are expected to run only logarithmically.) The tentative conclusion of our analysis will be 
that the position of the pole of a particle with mass comparable or larger than Planck's 
mass can be shifted to exponentially large energies. The sole remnant at lower energy 
would be the graviton. (In the generalized gauge theory of gravity alluded to above, one 
would also have an unbroken GUT sector). 

This paper will be organized as follows: in Section 2 we will describe in detail the gauge 
theory of gravity that we are going to consider. In Section 3 we discuss the linearization of 
the theory around flat space and define the effective potential for the dilaton. In section 4 
we define the average effective potential. In section 5 we write the renormalization group 
equations for the average effective potential and study their solutions. In Section 6 we 
draw our conclusions. 
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2. Lagrangians 

We consider a theory of gravity with independent metric g^^,, connection / a'^i/ ^ind dilaton 
p. The connection is required to be metric: 

dxg^iv - ^Qtu - rx'^i^QfiT = , (2.1) 

but can have nonvanishing torsion Tx^^,^ = Fx^u — r^^^x- The curvature of the connection 
will be denoted R^i^^a = d^^Fj^P cr — d^r^Pfj + Ff^Prrv^ a — r^^Prr^'^ a- It is antisymmetric 
in |U, z/, and using (2.1) one can show that R^vap = QapRiJ.y'' p is also antisymmetric in a, 
13. It has no other symmetry property. 

The most general diffeomorphism invariant action which is at most quadratic in the 
derivatives of the fields is 

S{g, r,p) = J d^x^/d^[^bo^^p^''p + 2b^pd^pTi' + V{p) + gop^R 

+ giR/.^paR""'"' + g2R^.paR''""' + gzR,..paR'^'"' + qaR^.R'"'' + g^R,..R'''' + g^^R^ 

(2.2) 

where R^j^ = Rp^^'v, R — g^'^R^if and T^, = T^.^^ ^. Indices are raised and lowered with 
g. The couplings of the dilaton can be understood as due to the quantum mechanical 
breaking of Weyl symmetry in a classically Weyl invariant theory. This point is discussed 
in Appendix A, but is not necessary for what follows. 

There are some special choices of coefficients that should be pointed out since they 
will be useful later. The first of these can be expressed by the identity 
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_ 1 /• 4 1 a/3,. ,.p. _ ^^-^^ 

" 128^ J d'x-j==e^^^'e^'''P''R^f,^,R^spa = X , 

^^1234 _ "1^^ which is an integral representation of the Euler invariant x- With this choice 
of action, the theory is topological. Another useful identity is 

R = R{g) + ^T^.pT^"" + ^T^vpT^^'" - T^Tp + 2V^T^, (2.4) 

where R{g) denotes the Ricci scalar of the metric g^jy and is the covariant derivative 
with respect to the Levi-Civita connection. Using this formula one can replace R by R{g) 
in (2.2) at the cost of redefining the constants aj and b^. In particular, if we choose the 
only nonzero coefficients to be bo, go, ai = —jgo, 0^2 = ^^go, 03 = ^0 and bi = 2go, (2.2) 
becomes the action for a nonminimally coupled scalar field: 



J d^x \/ det g 



Uod^pd''p + goR{g)p' 



(2.5) 
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For p = constant, this is Einstein's action. 

Finally, we mention two alternative ways of writing the torsion squared and curvature 
squared terms. In discussing the linearized theory it is useful to write them in the more 
compact form 

^M-/3p-75^^^^^ + p'i?'^""^^" V.T.^a , (2.6) 

where 

Qt^uadpa^s ^ SYMMj^i gf^pg^'^g^^gl^^ + g^ Q^'f g^" g""" Q^^^ + ^3 g'^'^ g^"" g'^" g^" 

+ ^"V(^4<?^V + ^75/V + ^76rt"')} , (2.7a) 
jj„auppa ^ SYMM{g^P{ai g'^'^g'"' + 03 i?""/^ + 03 £?"^/^)} . (2.76) 

The prefix "SYMM" indicates that one has to take the proper combinations so that G 
is antisymmetric in (/U, z/), (p, cr), (a,/3), (7, ^) and symmetric under the simultaneous 
interchange of fiuafS and pa'jS and similarly H is antisymmetric in {pjv), {p,cr) and 
symmetric under the simultaneous interchange of pau and p(3a. 

The other way of writing the action is based on the decomposition of the torsion 
and curvature tensors into irreducible parts with respect to the Lorentz group. This is 
discussed in Appendix B. 

3. Linearization 

In this section we set up the formalism for computing, at one loop, the Euclidean effective 
action Seff{g, F, p). We will not ultimately do the calculation (this was done for example in 
[11]), but what we describe here is a preliminary material for the definition of the average 
effective action, that will be given in the next section. We will restrict ourselves to a flat 
background 

gp.u = Sfiu , rx^iy = , p = constant . (3.1) 

In this case we can write Seff{g, T, p) = J d'^xVesip), where V^fj is the effective potential 
for the dilaton. The first step in this calculation is to linearize the action around the 
background (3.1). We define the fiuctuation fields SFx^i,, ^g^j^u and 5p = cr as the deviations 
of r, g and p from their background values (3.1). In principle there is a total of 64 + 10 + 1 
fields. However, the condition (2.1) implies that 

dxSg^u - SFx^y - STxufj. = . (3.2) 

(From now on indices are raised and lowered with 5pjy.) These are 40 constraints, reducing 
the number of independent fields to 35. 

We define u>xij.v = ^Txi^j^] and Lp^^, = ^p5g^u- The rescaling of Sg^i, is convenient for 
dimensional reasons: with this definition all the fiuctuations have dimension of mass. This 
redefinition is legitimate here since we are assuming p to be constant; it is related to the 
choice of functional measure in the path integral. We will see later that it does not have 
any effect on the renormalization group equations. After having written out the linearized 
action in terms of SF, Sg and a, one can use (3.2) to eliminate the symmetric part of SF in 
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favor of SFxi^^^j^) = ^^a^'/ui/- At this point the hnearized Euchdean action is a quadratic 
form which can be written, after Fourier transforming (we use 5^ — > iq^) 

S^^\^,uj,a;p)= \ [ d\Y. ^A{q)-OyAB]-^B{-q) , (3.3) 

^ A,B 

where the indices A, B label the three types of fields $i = a;, $2 = <^ and = cr and the 
dots stand for contraction over the tensor indices. When written out explicitly in terms 
of the components of the fields, O is a 35 x 35 matrix. The components of this matrix are 
listed below: 

C)M;.a/^'= 8(:;/«/^^^9,9, + 8iy^«/^V + 2^op'<^?(<^aM<^'^- , (3.4a) 

0[^^]a/^' = -SipH-^/^'q^ + 2igop [d^^d'^q^ - S^^q^ + 5'^5^Pq^) , (3.46) 

C[^^]a/" = 8iy^,^^^-g,g, + ^(^CM^^^ - 25^5^) , (3.4c) 
O^^^f''' = 2ipih - 2go)5''^q' , (3.4d) 

p dp 

= boq' + . (3.4/) 

For the purpose of clarity we have omitted to indicate explicitly symmetrizations and 
antisymmetrizations on the r.h.s. (for example (3.46) should be symmetrized in a, p and 
antisymmetrized in 7, 5). When V = 0, this linearized action is invariant under linearized 
gauge transformations. Let x'^ = x'^ — v'^ be an infinitesimal coordinate transformation. 
The variations of the fields are 

SFx", = dxd^v^" , = d^v, + d,v^ , 6p = 0. (3.5) 

There follows that the fields 

1 i 
^\,xv = ^qxiQixVu - qvV^) , (p^v = -p {q^v^ + q^v^) , ex = (3.6) 

are null vectors for the operator O. This can be verified by explicit calculations. To make 
O invertible, we add to the linearized action the following gauge-fixing term: 

5gf = ^ y d'^xd^v^'^dP^p, . (3.7) 

The ghost contribution has the form 

S^i^osMd) =\j d\d^0^dd];d, , (3.8) 



where, d, d, are anticommuting ghost fields and 



(3.9) 



Apart from the overall factor p, which can be eliminated by a redefinition of the measure 
and is irrelevant, this operator is field- independent. It can be neglected in what follows. 

To compute the one-loop effective action one now needs to calculate the functional 
determinant of the operator O appearing in the previous formulas. The determinant of O 
on the 35 dimensional space spanned by the fields is very hard to compute as it stands. 
One can partially diagonalize these operators in blocks corresponding to spin and parity. 
This is because C is a Lorentz covariant wave operator and therefore does not mix fields 
with different spin and parity. There are two modes with spin-parity 2+, coming from co 
and (f>, one 2~ mode from u;, two 1"*" modes coming from iv, three 1~ modes, two from u: 
and one from (p, four 0+ of which one comes from uj, two from ip and one from a, and finally 
one 0~ mode from u>. One counts indeed 2x5-|-lx5-|-2x3-|-3x3-|-4xl-|-lxl = 35. 

The total linearized quadratic action, including the gauge-fixing, ghost and potential 
terms, can be rewritten as 

5(2) = ljd'<l ^A{-q) ■ a^^iJ^) P,j^{J^) . ^B{q) , (3.10) 

where P^^{J^) are spin projection operators [12,8,10] that we list in Appendix C, and 
afj^{J^) are coefficient matrices, representing the inverse propagators of each set of fields 
with definite spin and parity. For V = these matrices are given by 



a(l+) 
a(l-) 



Giq^+Bip^ -iV2Bi\q\p 
iV2Bi\q\p B2q^ 



a(2+) = 
a(2-) ^Gsg'+Sip^ 



\/2S7p2 i^B^\q\p 
Bgp'^ iBs\q\p 
iV2Br\q\p -iBs\q\p Bsq^ 



G^q^ + Bqp 
V2Brp^ 



G^q^ + Bap-" 
i^/2Bg\q\p 


i\/6Bii\q\p 



a(0+) = 

. i 

a(O-) = Geq^ + 5i3p2 



iy/2BQ\q\p —iy/6Bii\q\p 

Bioq^ VSBm^ 



6og' 



V3Si2g2 



(3.11a) 
(3.116) 

(3.11c) 
(3.11d) 

(3.11e) 
(3.11/) 
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where 





= 4gi 


Go 


= 4(7i 


G3 


= 4(71 


Ga 


= 4(71 


^5 


— ^yi 


Ge 


= 4t/i 


Bi 


= 2ai 


B2 


= 4ai 


Bz 


= 6ai 


Ba 


= 2ai 



5a2 
3a2 



B^ = Aai - 2a2 , 

Bq = 2ai + a2 + 203 - go , 

B7 = a3- go , 

Bg = 2ai + 02 + 03 , 

Bq = 2ai + 02 + 803 - 2go , 

Bio = 4ai + 2a2 + 603 , 

Bii = bi- 2go , 

B12 = 2bi , 

Bis = 8ai - 8a2 - 2go , 
Bo = bo . 



(3.12) 



We observe that if we did not redefine the fiuctuation of the metric and worked with Sg, 
the only effect on the coefficient matrices would be to multiply by p/2 the second row 
and column of a{2'^), the third row and column of a{l~) and the second and third row 
and column of a{0~^). This would change the determinants of these matrices by an overall 
power of p^, which, as we shall see, does not affect the renormalization group equations. 

There are a few checks that one can make on these matrices. First we observe that 
the matrices a(l~) and a(0+) are degenerate. The proportionality of the last two rows 
and columns of a{l~) and the vanishing of the third row and column of a{0'^) are direct 
consequences of the diffeomorphism invar iance. 

If we take the only nonzero coefficients to be (73, g^j = —Ag^, g^ = g^, corresponding 
to the action (2.3), the coefficient matrices are identically zero; this is because the cor- 
responding action is a topological invariant (actually zero, since we expand around flat 
space). If we take the only nonzero coefficients to be go, ai = —\go, 02 = ~^90i cbs = ^O; 
bo and bi = 2go, corresponding to the Lagrangian (2.5), the coefficient matrices reduce to 



a{2+) 



a(0+) 





-2goq\ 

■0 

Agoq^ 4v^(7oq2 



.0 4VSgoq^ Boq^ 



(3.12a) 



(3.126) 



Note that for ^fQ = 6o/12 the second and fourth rows and columns of 0(0"*") are proportional, 
and the matrix has rank one. This is because in this case the action is Weyl invariant (see 
Appendix A). If we freeze p =constant, the last row and column of 0(0"*") can be suppressed 
and we are left with the familiar coefficient matrices of Einstein's theory [12]. 
The contribution of the potential to the inverse propagators is 



a(2+) 



Pot 







(3.14a) 
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a(l )pot = 



"0 











(3.146) 



a(0+)pot 



'0 








^/^ dV 
P dp 





p^ 

1 dV 
P dp 





V3 dV 
P dp 
1 dV 
P dp 
d^V 
dp^ 



(3.14c) 



Note that these matrices do not have the degeneracies of (3.11d, e) or (3.126). This is 
because flat space (with p ^constant) is a solution of the fleld equations only if F = 0. 
Finally the contribution of the gauge fixing terms is 



a(l )gf = 



a(0+)GF 















la ^ 





a. ^ 

OJ 



(3.15a) 



(3.156) 



With these results, the usual one- loop effective action is equal to the sum over spin 
J and parity V of the logarithms of the determinants of the matrices a. These are poly- 
nomials in qf^, p^, y and its derivatives of dimension up to eight. Taking into account the 
multiplicity of these contributions, the one-loop effective potential is 

-e.(.) = iB-.i)/|f.."(^^). (-) 

We have normalized the effective action with the determinant of a free field, which is 
obtained by linearizing the action aroimd a fixed constant field po- It is natural to choose 
Po as the minimun of Kfr itself, in which case 14ff(po) = 0. 

Given the previous explicit form for the matrices a(J^), one can now compute the 
expression for K!ff(p), using standard renormalization techniques. As explained in the 
next section, we shall follow instead a different strategy: we shall derive the equation that 
describes the renormalization group flow of Veff . This allows a more accurate discussion of 
the scale-dependence of the parameters that characterize the effective potential. 

4. Average effective potential 

To study the renormalization group flow of the effective potential, we shall use ideas 
originally introduced by Wilson [13]. One begins from some action S^-^ which is supposed to 
describe accurately the physics at some momentum scale k\ . Physics at a lower momentum 
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scale k2 is then described by an effective action which is obtained by functionaUy 
integrating exp(— 5"^^) over all fluctuations of the fields with momenta between ki and 

The procedure is then iterated: the effective action at scale ks < /c2 is obtained by 
functionally integrating exp(— ^^j) from k2 to ks, and so on. The flow of Sk with k is 
the renormalization group flow. Each step of the integration should not cover too large a 
range of momenta. In this way one can accurately compute the effective action at some 
low momentum scale k, starting from a high momentum scale A. Note that this is not the 
same as performing the functional integral from A to A; in one step, because the action is 
updated at each step of the integration. In the following we will refer to this updating as 
"the renormalization group improvement" . 

Originally this idea was applied to spin fields on the lattice, but it was subsequently 
adapted to the continuum, where it was used to clarify and simplify the notion of renormal- 
izability [14], and was also applied to gauge theories [15]. The particular implementation of 
this idea that we shall use here has been discussed in [16,17]. One can start from the usual 
definition of the effective action S'eff, defined as the Legendre transform of W = — InZ, 
where Z is the partition function. The effective action has a perturbative expansion in 
Feynman diagrams, which correspond to integrals of certain functions of the momenta. By 
introducing some kind of infrared cutoff k in the propagators one obtains a new effective 
action Sk, which depends on the scale k. 

A way of implementing this idea in the path integral formalism is to add to the 
classical action S a term ASk that constrains the average of the field in spheres of 
radius ^ k~^ centered around the point x to be equal to a predetermined function (pix) 
(one can apply this discussion to the theory we are interested in by replacing the generic 
notation with F, g and p). The average of a field around a point x is defined by the 
convolution fk4^{x) = J y ^J dei gfk{x — y)(f){y), or, when g is flat, in Fourier space, by 
/fc0(9) = fk{q)(l>iq), where fk is the function 



with a and b are constants. This function interpolates between a gaussian, for 6 = 1, and 
a step function for 6 — > oo. 

In [10] we discussed this procedure in the context of a gauge theory of gravity where 
only the coefficients gi and ai were assumed to be nonzero. The specific choice of the 
term AS^ that we used there had the disadvantage that some of the propagators were not 
well defined. This was not important in [10] because these terms did not contribute to 
the quantities that we computed there. However, it could be a problem in more general 
calculations. For this reason we shall use here a simpler definition: we will assume that 
the term AS^ is defined in such a way that the only effect it has on the linearized action 




(4.1) 
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is to replace in the operators O the term with the function 



(and \q\hy ^/Pkiqj). The function Pk{q) approaches exponentially fast the function for 
> A;^, but goes to a constant {b = 1) or diverges {b > 1) for q'^ 0. Thus, replacing q^ 
by Pk in the propagators suppresses the modes with q'^ < A;^ and therefore has an effect 
similar to putting an IR cutoff at momentum k. This definition of scale-dependent effective 
action Sk is equivalent, at least at one-loop order, to the one given in [17]. 

It is quite clear that the considerations that were made in [17] for nonabelian gauge 
theories can be extended in a rather straightforward way to the case of gravity. There 
is one point, however, that requires some extra care: it is the definition of the absolute 
normalization of the potential. In the presence of gravity, the value of the potential at 
the minimum is interpreted as a cosmological constant. It affects the field equations, and 
therefore cannot be shifted arbitrarily. Furthermore, we see from (3.14) that it corresponds 
also to the mass of the graviton. In principle, the value of the potential at the minimum 
could depend on the scale, but one has to make sure that at least in the extreme IR limit 
this value be zero, to ensure that the graviton be massless. 

We take the following definition of the scale dependent effective potential for p: 

where ak{J^) are obtained from the matrices a{J^) given in (3.11,14,15) by replacing q^ 
with the function Pk{q) and pk is defined to be the minimum of Vfe. Note that with this 
definition Vk{pk) = for all k, so that the cosmological constant is actually zero at all 
scales. In this connection see also [18]. 

Finally, we observe that since the determinants det afc(J^) are functions of p^, it is 
consistent to assume that Vk is a function of p^. It will be convenient to define V' = -^r^- 

Then in (3.14c) we can write 1^ = 2^' and ^ = 4pV' + 2V' . 
5. Renormalization group flow 

The average effective potential obeys a renormalization group equation that is obtained 
from (4.3) by taking its derivative with respect to k. It has the form 

where the functional F comes from the r.h.s. of (4.3): 

F(Vy,V") = \Y^(2J^,) I ^^_L_4(deta.(j'')) . (5.2) 



11 



with 

After taking these derivatives, we substitute in F the classical potential V with the effective 
potential 14: this is the "renormalization group improvement". This substitution gives 
(5.1), a differential equation for the function 14 (p). Notice that thanks to the behaviour 
of Pk and its ^-derivative in (5.3), the integral in (5.2) is actually dominated by a finite 
range of momenta and docs not need an ultraviolet regularization. 

Although derived in the context of a one loop approximation, this renormalization 
group improved equation has a validity that goes beyond one loop [19,13,14]. Clearly one 
cannot follow the evolution of the whole function T4 ) so some other kind of approximation 
becomes necessary. In the following we shall study only the first few terms of the Taylor 
expansion of Vk- As explained in the introduction, we are interested in the spontaneously 
broken phase, with a nonzero VEV of p aX k — 0. Thus, we parametrize T4 by the position 
of its minimum, p^, and the quartic coupling at the minimum Xj.: 

Vlipk) = , Afe = Vl'ipk) . (5.4) 

In terms of these parameters the potential reads 

Vkip) = IXkip' - pI)\ (5.5) 

which is the Taylor expansion of T4 cis function of around its minimum. 

The equations governing the evolution in k of p^ and can be obtained by differen- 
tiating the definitions (5.4). This leads to the following set of coupled partial differential 
equations: 

k^=l{k)k\ (5.66) 



dk 

k^=p{k), (5.6c) 



with 



k^n) (pk) = ^/ dxxn^{Pk,pl)k^ , (5.76) 

^^^^ = {^^k^^) ^P""^ = 3^1 dxxnpiPk,pl)k^ . (5.7c) 

In the definition of (3{k) we are neglecting a term containing V^"{pk), which takes into 
account the k dependence of the point of definition of and is a signal of the presence 
of operators of higher dimension in the potential ("irrelevant" operators). Consistent with 
our approximation, we will neglect the effect of these terms. The functions (3 and 7 have 
the general form shown in the r.h.s. of (5.7), where x = and are rational functions 
of dimension k~^ . These functions can be computed explicitly from the representation 
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(4.3) of Vfe, by taking derivatives with respect to and k, and using (5.3). The general 
expressions are comphcated and not particularly illuminating, so we will not give them 
here. It is not possible to find a solution of the system of p.d.e.'s (5.6-7) in closed form. 
However, analytic solutions can be obtained in some asymptotic limit. 

Let us assume first that pi is small with respect to k. One would expect this to 
describe the behaviour of the theory in the regime when k is large compared to the Planck 
mass. In simpler systems, this approximation indeed reproduces the results of perturbative 
calculations of beta functions at momenta much larger than the characteristic mass of the 
theory. Because of the factor k^^, the integrals are dominated by the region x ~ /c^. 
In this region, is itself of order fc^, so that is small with respect to Pk- One can 
therefore expand the functions TZ in powers of p|. In the function TZj the dominant term 
is a constant (independent of pk) but the function 7^^ has a pole for p^ 0, coming from 
the contribution of the spins 2+, 1~ and 0+. The beta functions reduce to the following: 



7(/^) = 70 + O(g 



Pi 



p / 2 

/3(A;)=/3_i^+/3o + o(^|| 



(5. 



where 



7o = 



327r2 Afc 

6Afc(-Bio 



B1B2 - 2B^ _|_ _|_ 3 -^3-^5 
G1B2 G2 G^B^ 



GaBs 



+ 



2B 



12; 



BaB 



O-DlO 



p-1 



1 



327r2 
1 



10 

B~2 



B9 

H ^ - 

Gs 

Bo 



327r2 



BqBiq — 35^2 
?2^2 



2BqBq — 12BgBiiBi2 + QBioBfi 
G5(BoBio — 35^2) 

+ 13q; Afc/_2 • 



+ 



B 



13 



Ge 



1-1 



2\2 



10 



{B^B2-2Blf + A\uBl ^ ^Bf ^ (BsB^-2Bl) 



GlBl 



+ 10^+6- 



GlBl 



, ^ [BqBs - 2^2)2 24(25oSio + 3S2o - I2S10S12 + 7^22) - 16a5io -2 



n2 d2 

432(Sio5ii - SgSia)^ 



{BqBio — 3523)2 
4{BoBg — 3BiiBi2 + I2S10S11 — 1259^12)2 



^5(^0510- 3522)2 

2{2BqBq — BqBqBio + 6S10-B11 — 12BqBiiBi2 + 3BgBi2)^ 



Gl{BoBio — 3512)2 
The positive dimensionless constants are defined by 

dPk 



+ 2 



Biz 
Gl 



1-3 = 

1-2 = k- 



o t 

dxxP^. /c- , 
dk 



1 



dxx Pi, k —7—; 
" ok 



(2a) 1/6 



r(l + 1/6) 



dxxk 



dPk 
dk 



(2a)3/& 



r(l + 3/6)C(3/6). 



(5.9a) 
(5.96) 



Afe (5.9c) 



(5.10) 



Infrared convergence of Iq requires that 6 < 3. 
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Note that in this approximation the dominant term in P{k) is the one coming from 
the pole, unhke other known examples where the dominant term is the constant Pq [16]. 
This peculiarity can be traced to the presence of the undifferentiated potential V in the 
inverse propagators, which is characteristic of gravity (cfr. (3.14)). The situation would 
be different if we allowed a nonzero cosmological constant. 

Treating 70, and /?o as constants and neglecting integration constants one finds 
that p1 = |7o/s^ and Afc = 2^^1n/c^. This is the behaviour that one would expect on 
dimensional grounds. In this calculation one neglects the running of the coupling constants. 
This is a reasonable approximation if one considers the behaviour of the theory over a range 
of momenta which is not too large. One could take into account the running of Afe, for 
which the evolution equation is known, but the result would not be very significant: the 
other coupling constants are also expected to run logarithmically, and their (unknown) 
contribution could easily overwhelm the one coming from Afc. 

Assuming that all the couplings appearing in 70 run logarithmically, the true solution 
for p^. would deviate from the one given above by a sublogarithmic correction. The validity 
of the approximation p| <S k'^ hinges on the sign of these corrections. If p| grows slower 
than /c^, the approximation could be justified. This was the case in the calculation we did 
in [10]. However, it may not be generally true. 

Given that, in general, the validity of the approximation pf, <C is questionable, 
it would be desirable to say something on the large-A; behaviour without making this 
assumption. One general conclusion that can be drawn with reasonable confidence is that 
p| will be proportional to fc^, up to (at most) logarithmic corrections. To see this consider 
again the general form (5.7) of the functions 7(/c) and P{k). As mentioned before, the 
x-integration is cut off exponentially for x > k'^, and as a power for x < k^, so if we are 
only interested in the dominant behaviour of the integral, we can replace Pf- by A;^. Assume 
further that pi = ck'^ for some constant c. Then 7?.(fc^, c/c^) is a constant that can be taken 
out of the x-integration and the functions 'j{k) and f3{k) become simply constants. The 
equation for p| becomes an algebraic equation that implicitly determines the constant c, 
so the assumption = c/c^ is justified a posteriori. 



Let us now consider the opposite limit: k <S 



pI 



This is the limit A; — > 0, when 



po 7^ 0. In this case we retain in each sum the term with the highest power of p^. The 
contributions of the different spins are not all of the same order in /c^/p^. Keeping only 
the leading terms, the beta functions reduce to 

k'^ / 

7(/^) = 72^ + o(- ; Pik) 



Pk 



4 



+ 



Pi 



(5.11) 



where 
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8XkBQ{2B9 — Bio)' 



B 



13 



14 



10 _ 5 
B2 — 2Bi 2Bg — BiQ 



+ 15a 



(5.126) 



In this regime the running of Afe and is damped by powers of /c^/p^ and stops for /c — ^ 0. 
The solutions for smaU k are 



Pk 



pI 



^k — ^0 



1 + 



1 + 



72 


k^- 


6 


pI. 


Pi 




2 


Pi. 



(5.13a) 
(5.136) 



For generic values of the parameters Qi, ai, bi, all modes except for the graviton are 
massive, with masses of the order of Planck's mass. One would expect that these modes 
can be neglected when describing the physics below the Planck scale. Thus, the running of 
the dilaton potential at low energies should be derivable entirely from the action (2.5). This 
can be easily checked using the coefficient matrices given in (3.12). The renormalization 
group equations become for small k 



k 



k 



dpl_ 3 12^0 
dk 

dXk 



327r2 
15 

dk ^32^2' 



8AI 



So, k 
— Jo 



6 



Pk 



^k{a - — 



1 

S90 



I ^ 



(5.14a) 
(5.146) 



These have again the solutions (5.16). Note that the coefficient 72 vanishes when ^0 = 
60/12, in which case the action (2.5) is conformally invariant. (See Appendix A). 



6. Concluding remarks 

One attractive feature of the theory we have considered here is that one can easily describe 
the transition to Einstein's theory at sufficiently low energies. Indeed, for generic values 
of the coupling constants Qi and a^, all components of the connection F^'^ p are massive, 
with squared masses of the form p^B/G^ where B and G are appropriate combinations 
of the coupling constants and gi. As mentioned in the introduction, this is due to the 
occurrence of a Higgs phenomenon. The dilaton becomes massive too, with a mass squared 
equal to AoPo- energies lower than the VEV po, all these particles decouple, leaving 
the graviton as the only remnant. Its effective dynamics is given by the second term in 
(2.5). If we assume that the Einstein term of the low energy world comes entirely from 
this source, we see that po has to be of the order of the Planck mass, mp. 

The results of the preceding section can now be summarized as follows: at scales 
k much larger than mp, the coupling constant runs logarithmically and the VEV p^ 
runs quadratically, as one would expect on dimensional grounds. These results are the 
reflection of the logarithmic and quadratic divergences that one would encounter if one 
tried to remove the ultraviolet cutoff from the theory. On the other hand for scales k 
below the Planck mass the running is suppressed, and both A^ and pk tend to constants 
for /c ^ 0. 
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When we said that there are massive particles in the theory we used in the mass 
formula the VEV of p at the scale = 0. This is the naive procedure that one would 
follow, at tree level, with a classical potential. However, the running of the VEV of the 
dilaton affects the position of the pole of the propagators. Let us see this first in the case 
of the dilaton itself. From the parametrization (5.5) of the potential T4, the running mass 
m?{k) = ^y-^l is equal to A^p?. The physical mass of the particle, i.e. the position of 

' Pk 

the pole of the propagator, is defined implicitly by the condition 

™phys = ™'WL2=^2 • (6.1) 

phys 

This is equivalent to the statement, which was made in [10], that in a propagator of the 
form {q^ — m^(/c))~^ the pole is to be found by evaluating the running mass at k = \q\. 
Equation (6.1) is solved graphically, by finding the intersection of the plots of the l.h.s. 
y = 'Tiphys and the r.h.s. y = XkPk\ j^2-^2 , as functions of the independent variable 

phya 

mpjjyg. The l.h.s. is represented by the straight line at 45*^. The r.h.s. is the plot of 
the running of pk that we have computed in the previous section, multiplied by a factor 
Aq, up to logarithmic factors that we neglect for the moment. As we have seen, it starts 
fiat at AoPo and grows linearly for large Wphyg- If Aq <^ 1, the intersection occurs in the 
region where the r.h.s. is constant, and therefore the mass squared is AoPo ^ ^p- 
the other hand if Aq is of the order or bigger than one, the intersection is shifted to higher 
energies. Exactly where it occurs depends crucially on the logarithmic corrections. If the 
leading logarithmic factors in the r.h.s. appear with a negative power, Wpj^yg is of the order 
exp(Ao)po- If they appear with a positive power, there may be no intersection at all. In this 
case, the dilaton would disappear completely from the spectrum. Note that an anomalous 
dimension of p would give a power correction to the running of the mass, so it would be 
even more important than the logarithmic corrections in the previous considerations. 

This discussion can be repeated for the connection Fx^^, which potentially carries 
dangerous ghost or tachyon states. The running of pk could eliminate these states from 
the spectrum. That such a mechanism could exist was suggested in [20], but no concrete 
support for this idea had been given until now. In order to draw some definite conclusion 
one would have to compute the beta functions for the couplings gi and for large k, find 
the ultraviolet fixed point, if there is any, and evaluate (5.9) at that point. 

The same discussion can be repeated also for any other matter field. In this theory the 
matter is supposed to be coupled to gravity in such a way that the masses of all particles 
are proportional to the VEV of the dilaton. For example, in the case of a scalar field 0, 
the action would be 

S{cf>,g,p)=^ J d^xy/d&rg[g^^d^4>dA + hp^<t>^] , (6.2) 

where /i is a dimensionless coupling constant. The reason why the known particles have 
masses much smaller than Planck's mass would be the smallness of the coupling constant h. 
For such particles the intersection of the curves given by the l.h.s. and r.h.s. of (6.1) occurs 
in the region where the VEV of the dilaton is constant, so the poles are exactly where one 
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would expect to find them. In GUT theories the expected masses are only a few orders 
of magnitude smaller than the Planck scale, so that the renormalization group corrections 
envisaged here may not be negligible. In the case of a scalar field, this is discussed in a 
separate publication [21]. 

Finally, we mention that the definition (4.3) of the scale dependent effective potential 
is not the only possibility. One could choose a different normalization of Vk{p) such that 
its minimum (the cosmological constant) is actually /c-dependent. These alternatives could 
be of relevance for example in cosmological problems. 
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Appendix A: Broken Weyl invariance 

Let us consider the behaviour of the action (2.2) under Weyl scahngs 



^M-^^V > p^n-V , rx^,^rx'', + si^n-^dxn , (ai) 

where is a general function of the position. There is also another possible way in which 
a Weyl transformation could act on the connection: it is defined by requiring that F trans- 
forms like the Christoffel symbols of g. This alternative transformation acts trivially on the 
torsion tensor, whereas [A.l) acts in a nontrivial way. We have chosen the transformation 
(^.1) because it can be generalized to a local GL{4) transformation. 

It is easy to see that the curvature tensor R^.^'^a is inert under (^.1), so the curvature 
squared terms are invariant (as in Yang-Mills theory). The remaining terms can be made 
invariant by choosing 

6o = 6(2ai + 02 + Sas) , 6i = 2ai + 02 + 803 , V{p) = ^Ap^ , {A.2) 

where A is an arbitrary constant. Note in particular that with the choices leading to (2.5) 
this gives go = bo/12. 

It is sometimes convenient to write the action in a different form. Define the combi- 
nations 

~9^.u=p'g^.. , rx''u = rx^v + Kp~^dxp , (A3) 

which by construction are inert under the transformation {A.l). We also define the cur- 
vature of f, R^uf^a = d^r^P^ - d^r^P^ + r^Prh^a - r^^rr^/a, and the torsion of T, 
fx^'u = rx^u - ru^x- The following relations hold: 

fx\ = Ta^ + ^^^p-'^AP - 5>-'a,p , iA.4b) 

dx9,,u - rx'^^.gr^ - rx'^j^g^r = P^ {^x9^^u - rx'' ^,grv - ^x''vg^,T) ■ {A.Ac) 

Note from {A.Ac) that if p 7^ 0, metricity of one connection guarantees the metricity of the 
other. 

With the relations (^-2), the action (2.3) can be written S{g., F, p) — S{g, F), where 



S{g,F)^ J d^x a/ det g 

+ giRt^upaR^^"'"' + 92R^^.paR'"''"' + 9^R^..paR'""''' + 9^R^..R'''' + 95R^..R''^ + ^6^'] 

where R^^, = Rp^P^, R = g^^R^u and Ti, = Tj^^ In this action indices are raised and 
lowered with g. This is the most general action quadratic in T and R, and it is invariant 
under the transformations {A.l) in a trivial way, since all quantities entering in this action 
are inert under those transformations. It is obtained from (2.2) by choosing the conformal 
gauge in which p is constant, and rescaling the fields by factors p. Conversely, (2.2) with 
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the relations {A.2) is obtained from (^.5) by using the definitions {A.3). Note that the 
field p, which we call the dilaton in this paper, may be called the conformal factor of the 
metric g if one started from the action (A.b). In the literature on conformal gravity the two 
actions (2.3)-(A.2) and {A.5) are said to be written in the Jordan and Einstein conformal 
frames, respectively. They are completely equivalent at the classical level. 

Let us now consider the quantization of the theory in the Weyl-invariant case. We 
note first of all that when (^.2) holds, Bn = Bg and B12 — -Bio- Thus at the linearized 
level the new gauge invariance manifests itself in the proportionality of the second and 
fourth rows and columns in a{0^). Therefore, one has to fix the gauge also for Weyl 
transformations. Assuming that this has been done, it is easy to see that there exists a 
quantization procedure that preserves Weyl invariance. One has to use the form (^.5) for 
the action, and define the functional measure by means of the metric g. In this case the 
effective action can be written again as a functional of g and F alone, and therefore is 
automatically invariant under (A.l) [22]. 

One is naturally inclined to preserve as much as possible the classical symmetries in 
the quantization process, so this choice of measure may seem to be the only sensible one. 
However, this is not the case. Other choices are possible and, from a certain point of view, 
may even be more natural: if we interpret the metric g as the one defining the geometry 
of spacetime, then it is natural to use g rather than g in the definition of the functional 
measure, and this leads to a quantum theory in which Weyl invariance is broken. 

In a concrete calculation, the definition of the measure refiects itself in the definition 
of the cutoff. We have shown in [10,11], that if one starts from a Weyl invariant theory 
and defines the UV momentum cutoff as g'^'^q^Qi, < (rather than g^^Qf^qv < A^), then 
the effective action Sef[{g, F, p) will not have Weyl invariance anymore. 

In this paper we have only discussed the renormalization group fiow of the effective 
action, which does not necessitate the explicit introduction of an UV cutoff. Still, Weyl 
invariance cannot be maintained. This can be seen as follows. Suppose we study the small 
fluctuations of the gravitational field around flat space at some energy scale ki , and suppose 
that these are well described by the Weyl invariant action (2.3)-(A.2). Now suppose we 
want to know the effective action at some lower energy scale /c2- As discussed in section 
4, this is given by a functional integral over all fiuctuations of the fields which lie in the 
momentum shell between ki and k2. But how is this momentum shell defined? Since we 
are postulating that the geometry of spacetime is given by the dimensionless metric g, the 
shell is defined by k2 < g^^q^qv < ki (for simplicity we are assuming here a sharp cutoff, 
but this is by no means essential). One sees that the definition of the shell introduces in 
the definition of the effective action a dependence on g which is not compensated by a 
dependence on p. Unavoidably, will not be a function of the combinations (2.4) alone. 
It will be a genuine functional of g, F and p, and will not be invariant under infinitesimal 
Weyl transformations. 

To be more specific, suppose that we want to compute only the scale-dependent effec- 
tive potential Vk{p), which is given by 5"^ = / d^x-^det g Vfc, for a field of the form (3.1). 
As we have mentioned above, Weyl invariance requires the potential to be purely quartic. 
If we assume that V^j is purely quartic, then unavoidably Vk^ will not be anymore, because 
the integration procedure breaks scale invariance. It is instructive to see this in detail in an 
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explicit calculation. Consider a single wave operator of the form O = + cp^ (in principle 
(3.3) can be written as a sum of such terms), and suppose the momentum shell is defined 
by sharp UV and IR cutoffs. Then we have 



1 
2 



647r2 



fc| In 



kl + cp^ 



(2^ 



In 



q^ + cp^ 
g2 + cpl^ 



+ c p m ' 



kl + cp2 
kl + cp^^ 



which is obviously no longer purely quartic. It is still a function of p^ only. 
Let us approximate the potential T4 by a quartic polynomial of the form 



(A.6) 



14 (P) = -jTulp^ + Xkp'^ . 



{A.7) 



The constants m| and Afc are to be thought of as the coeflBcients of the Taylor expansion 
of Vk around the origin: 



2 d'Vk 



dp' 



Afc 



1 d^Vk 



24 dp^ 



Then we have 



Ami = ™fci 



c 
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fe2 



327r2 

^2 



^Afc = Afc^ - Afc, 



{kl kl) 
kl 



(^.9a) 
(AM) 



The most important conclusion of this discussion is that the potential cannot be 
consistently assumed to be purely quartic at all scales: if is zero at some scale, it 
will be nonzero as soon as one begins to integrate. Weyl invariance is broken and a mass 
term is generated. There follows that if we want to study the renormalization group flow 
of the effective potential taking into account what we called the "renormalization group 
improvement" , we have to assume from the outset that the potential is not purely quartic. 
Of course, if we were to study the whole effective action rather than just the effective 
potential we would flnd many more terms that are not present in the starting action. For 
example, instead of the factors p^ in front of the torsion terms there will now be general 
functions of p. In principle, this will also have an effect on the running of the potential, 
but we neglect this effect. 
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Appendix B: Lorentz decomposition 

In the main text we have taken as independent couphng constants the coefficients of all 
possible contractions of two curvature or torsion tensors. However, these contractions do 
not carry any special geometrical significance. One may expect that the final results would 
look simpler if they were expressed in terms of another set of parameters, which are related 
to geometrically significant quantities. Let us consider the irreducible parts of the torsion 
and curvature tensors with respect to the Lorentz group [23] : 
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{9pp9ua - QpaQup) R , (S.3a) 
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22 ^9pp9va 9pa9iyp) 



(S.36) 
(B.3c) 

(B.3d) 
(S.3e) 
(S-3/) 



where R^^'^^i, = \ {R^y + R^p) - \9piyF and R^^\v = | (Rpi' - Rup) are the symmetric 
traceless and antisymmetric parts of the Ricci tensor. These decompositions are orthogonal 
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with respect to the inner products {A,B) = J d'^x^^det gg^'^g'^^g^'^ ■ • ■ A^j^\___B 
quantity (2.6) can then be rewritten in the form 
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where and gi is a new set of couphng constants, related to the couphng constants a^, gi 
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(S.5) 

Most formulae are more compactly written in terms of the new coupling constants. This is 
evident from the following table, which gives the relation between the coefficients Gi and 
Bi appearing in (3.11) and the parameters appearing in the Lagrangian: 
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Appendix C: Spin-projector operators 

For completeness, we list in this appendix the explicit expressions of the spin-projector 
operators P^^{J^) that have been used to rewrite the action (3.3) into the form (3.10). 
There are in fact some differences with respect to those in [10], due to the different spin- 
parity content of the fields appearing here. 

For fixed spin J and parity V, these operators are labelled by the indices A, B that 
identify the fields a;, <^ and cr, and by z, that identify isomorphic Lorentz representations 
occurring more than once. For example for spin-parity 2"*", i = 1,2; for 1~, i = 1,2, 3 etc. 
The operators P^'^{J^) project out of a field a given irreducible representation, while the 
intertwiners P^^{J^) (with i ^ j) give isomorphisms between the different representations 
occurring more than once. (Note that the indices A, B in these projectors are redundant 
since z, j already label the representations. It is nevertheless convenient to keep them in 
order to remember by what field a certain representation is carried, e.g. for = 1~ the 
representations i = 1 and i = 2 are carried by co and z = 3 is carried by <f). The operators 
P^^{J^) are orthonormal and complete: 



P^ 



AB 



(J^) 



J,'P,A,i 



It is useful to introduce the following notations: 



(C.l) 



obeying the relations: 



(C.2) 



fj, " 



T P 



In terms of q^, L'j^, and T^, one finds: 

[^iT(2+)].p.'^^^ [pr2^{2+)i,r^ 

[^^27 (2+)],. [^^27(2+)],."^ 



[P(2+)] = 



rP 



(C.3) 



[p-(2+)] 



a/37 



rpa 



= T? + Tp Lfj T-^ - - r,[, Li] T^]« , 

x/2, 



K^(2+)],,,"^ = ^T(«T[^;g.] - ^T^/'T.ipg.] 



[p--(2-)] "/'^ = \vt1^t2\ + \tPt^}k^-t^[pT'^^t^^- , 
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JTpcr 
] 



[pr(i+)] 
[PiT(i+)] 



q:/37 T-ia rp\P 7-7] T^[7 r/^l 

rpa - J-[p ^a] " ^[p-'a] ' 



a/37 



rpa 



■ a r7^[/3 ^7^7] 



[^'iT(i-)].p."^^ 
[P2T(i")] 



a/37 



It pa 



[^2T(i-)].p."^^ 
[^'3T(i-)]p.'^^^ [^^3T(i-)]p."^^ 



iTpa 



[p,7(i 

[^27(1 



a/j - 

rpo" 

a/3 

rpcr 



JpCT 



[p,7(i- 

[^2T(l" 
[^^2^(1" 

[p,7(i- 

[P37(l- 



1 a/37 7^ rpll rp(3]a 

J rpa '''[P cr] ' 

].p."'" = ^^r[P^]^''"' 
]^^^a/37 = 2L,"L|^T;;, 

] 



Tpa 



af3j 



I pa 
pa 
pa 



«/37 = 2r4^T;j , 



a/3 



(p 
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[P(0+)] 



[^4^(0+)]''^^ [P47(0+)]"^ 



irpcr 



[^27(0+)],."'^ 
[^^37(0+)],."^ 



ir pa 



[^27(0+)],. 
[^^37(0+)],. 
[^Y/(o+)] 



[Pi"i"^(0+) 

K70+) 
K70+) 
[PiT(o+) 

[pr(o+) 
K70+) 
K70+) 

[P37(o+) 
K70+) 

[P37(o+) 
[Pir(o+) 

K70+): 



a/3 



2 

rpa V 3 

x/2 



L''^ Trip ^a] 



pa 



J- pa J- y ) 



1 



a/3 rp rpa/3 

pa - Z ^'^ ' 

a/3 rp r a/3 

- J-T 



pa 



pa 



a/3 r rya/J 

nrr ^ •> 



pa 



a 13^ 



= L 



pa 1 



;[7 y/^la 



_ J_ ^a/3 



1 9 

^p'^^^^'n~^l Ck/37 T-ia T-i[/3 T-iT] rp[~^ rp(}\ rpot — rp[a rpP rpj] 

hrpa ~ 2 ^ ''^ 3 ^[p ^c^] - ^[r ^a] ' 
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